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When electrons are subject to a large external magnetic field, the conven- 



tional charge quantum Hall eff'ect |2l| dictates that an electronic excitation 
gap is generated in the sample bulk, but metallic conduction is permitted at 
the boundary. Recent theoretical models suggest that certain bulk insulators 
with large spin-orbit interactions may also naturally support conducting topo- 



logical boundary states in the extreme quantum limit [3|, , which opens 

up the possibility for studying unusual quantum Hall-like phenomena in zero 
external magnetic fields J 6|. Bulk Bii ^iSb^; single crystals are predicted to 



be prime candidates j?, Isl] for one such unusual Hall phase of matter known 
as the topological insulator 0, [lol . The hallmark of a topological insu- 

lator is the existence of metallic surface states that are higher dimensional 
analogues of the edg e states that characterize a quantum spin Hall insulator 
0, Q, Q, S, 0, y, y, 10, H, Q, 3]. In addition to its interesting boundary states, 
;he bulk of Bii_T^Sb,,: is predicted to exhibit three-dimensional Dirac particles 

findings of two-dimensional graphene and charge quantum Hall frac- 

tionalization observed in pure bismuth (211]. However, despite numerous trans- 



171, another topic of heightened current interest following the new 



port and magnetic measurements on the Bii ^Sb^ family since the 1960s 17|, 
no direct evidence of either topological quantum Hall-like states or bulk Dirac 
particles has ever been found. Here, using incident-photon-energy-modulated 
angle-resolved photoemission spectroscopy (IPEM-ARPES), we report the di- 
rect observation of massive Dirac particles in the bulk of Bio.gSbo.i, locate the 
Kramers' points at the sample's boundary and provide a comprehensive map- 
ping of the topological Dirac insulator's gapless surface modes. These findings 
taken together suggest that the observed surface state on the boundary of the 
Dulk insu 



ator is a realization of the much sought exotic "topological metal" 
111 ]. They also suggest that this material has potential application in 
developing next-generation quantum computing devices that may incorporate 
"light-like" bulk carriers and topologically protected spin-textured edge-surface 



currents. 
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Bismuth is a semimetal with strong spin-orbit interactions. Its band structure is beheved 
to feature an indirect negative gap between the valence band maximum at the T point 



of the bu 
points 



13, 



k Brillouin zone (BZ) and the conduction band minima at three equivalent L 



22j (here we generally refer to these as a single point, L). The valence and 



conduction bands at L are derived from antisymmetric (L^) and symmetric (L^) p-tjpe 
orbitals, respectively, and the effective low-energy Hamiltonian at this point is described 
by the (3+l)-dimensional relativistic Dirac equation Q, Q, [l^. The resulting dispersion 
relation, E(k) = -^\J {v -k) + v ■ is highly linear owing to the combination of an 

unusually large band velocity v and a small gap A (such that |A/|w'|| ?s 5 x 10~^A~^) and 
has been used to explain various peculiar properties of bismuth 14, Q, Q. Substituting 
bismuth with antimony is believed to change the critical energies of the band structure as 
follows (see Fig.l). At an Sb concentration of a; ~ 4%, the gap A between and closes 
and a massless three-dimensional (3D) Dirac point is realized. As x is further increased this 
gap re-opens with inverted symmetry ordering, which leads to a change in sign of A at each 
of the three equivalent L points in the BZ. For concentrations greater than x ~ 7% there is 
no overlap between the valence band at T and the conduction band at L, and the material 
becomes an inverted-band insulator. Once the band at T drops below the valence band at 
L, at X ~ 8%, the system evolves into a direct-gap insulator whose low energy physics is 
dominated by the spin-orbit coupled Dirac particles at L 

Recently, semiconductors with inverted band gaps have been proposed to manifest the 
two-dimensional (2D) quantum spin Hall phase, which is predicted to be characterized by 



ill pnast 

1m h, 



the presence of metallic ID edge states 3|, U, l5|, ll2|] . Although a band- inversion mechanism 



and edge states have been invoked to interpret the transport results in 2D mercury telluride 
semiconductor quantum wells no ID edge states are directly imaged, so the.r topological 

character is unknown. Recent theoretical treatments P, |8| have focused on the strongly spin- 
orbit coupled, band-inverted Bii_^.Sb^. series as a possible 3D bulk realization of the quantum 
spin Hallphase in which the ID edge states are expected to take the form of 2D surface 
states [7,y, 9] that may be directly imaged and spectroscopically studied, making it feasible 
to identify their topological order parameter character. Most importantly, the 3D phase is 
a new phase of matter in terms of its topological distinctions jlo| . 

High-momentum-resolution angle-resolved photoemission spectroscopy performed with 
varying incident photon energy (IPEM-ARPES) allows for measurement of electronic band 
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dispersion along various momentum space (fc-space) trajectories in the 3D bulk BZ. ARPES 
spectra taken along two orthogonal cuts through the L point of the bulk BZ of Bio.gSbo.i 
are shown in Figs la and c. A A-shaped dispersion whose tip lies less than 50 meV below 
the Fermi energy (Ep) can be seen along both directions. Additional features originating 
from surface states that do not disperse with incident photon energy are also seen. Owing 
to the finite intensity between the bulk and surface states, the exact binding energy {Eb) 
where the tip of the A-shaped band dispersion lies is unresolved. The linearity of the bulk 
A-shaped bands is observed by locating the peak positions at higher Eb in the momentum 
distribution curves (MDCs), and the energy at which these peaks merge is obtained by 
extrapolating linear fits to the MDCs. Therefore 50 meV represents a lower bound on the 
energy gap A between and L^. The magnitude of the extracted band velocities along the 

and ky directions are 7.9 ± 0.5 x 10^ ms~^ and 10.0 ± 0.5 x 10^ ms~^, respectively, which 
are similar to the tight binding values 7.6 x 10^ ms~^ and 9.1 x 10^ ms~^ calculated for the 
La band of bismuth j22|. Our data are consistent with the extremely small effective mass 
of 0.002me (where rrie is the electron mass) observed in magneto-reflection measurements 
on samples with x = 11% |23]. The Dirac point in graphene, co-incidentally, has a band 
velocity {\vf\ ~ 10^ ms~^) [18] comparable to what we observe for Bio.gSbo.i, but its spin- 
orbit coupling is several orders of magnitude weaker jsj], and the only known method of 
inducing a gap in the Dirac spectrum of graphene is by coupling to an external chemical 
substrate [20|]. The Bii_a.Sba; series thus provides a rare opportunity to study relativistic 
Dirac Hamiltonian physics in a 3D condensed matter system where the intrinsic (rest) mass 
gap can be easily tuned. 

Studying the band dispersion perpendicular to the sample surface provides a way to 
differentiate bulk states from surface states in a 3D material. To visualize the near-Ep 
dispersion along the 3D L-X cut (X is a point that is displaced from L by a fc^ distance of 
Sn/c, where c is the lattice constant), in Fig. 2a we plot energy distribution curves (EDCs), 
taken such that electrons at Ep have fixed in-plane momentum {kx, ky) = (L^, Ly) = (0.8 
A^^, 0.0 A^^), as a function of photon energy {hu). There are three prominent features in 
the EDCs: a non-dispersing, k^ independent, peak centered just below Ep at about —0.02 
eV; a broad non-dispersing hump centered near —0.3 eV; and a strongly dispersing hump 
that coincides with the latter near hu = 29 eV. To understand which bands these features 
originate from, we show ARPES intensity maps along an in-plane cut KMK (parallel to 
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the ky direction) taken using hu values of 22 eV, 29 eV and 35 eV, which correspond to 
approximate kz values of L^— 0.3 A^^, L^, and + 0.3 respectively (Fig. 2b). At hu 
= 29 eV, the low energy ARPES spectral weight reveals a clear A-shaped band close to Ep. 
As the photon energy is either increased or decreased from 29 eV, this intensity shifts to 
higher binding energies as the spectral weight evolves from the A-shaped into a U-shaped 
band. Therefore the dispersive peak in Fig. 2a comes from the bulk valence band, and for 
hu = 29 eV the high symmetry point L = (0.8, 0, 2.9) appears in the third bulk BZ. In the 
maps of Fig.2b with respective hu values of 22 eV and 35 eV, overall weak features near 
Ep that vary in intensity remain even as the bulk valence band moves far below Ep- The 
survival of these weak features over a large photon energy range (17 to 55 eV) supports their 
surface origin. The non- dispersing feature centered near —0.3 eV in Fig. 2a comes from the 
higher binding energy (valence band) part of the full spectrum of surface states, and the 
weak non-dispersing peak at —0.02 eV reflects the low energy part of the surface states that 
cross Ep away from the M point and forms the surface Fermi surface (Fig. 2c). 

Having established the existence of an energy gap in the bulk state of Bio.gSbo.i (Figs 
1 and 2) and observed linearly dispersive bulk bands uniquely consistent with strong spin- 
orbit coupling model calculations 



16 



22| (see Supplementary Information for full 



comparison with theoretical calculation), we now discuss the topological character of its 
surface states, which are found to be gapless (Fig. 2c). In general, the states at the surface 
of spin-orbit coupled compounds are allowed to be spin split owing to the loss of space 
inversion symmetry [E{k, f) = E{—k, f)]. However, as required by Kramers' theorem, this 
splitting must go to zero at the four time reversal invariant momenta (TRIM) in the 2D 
surface BZ. As discussed in [?, 9|, along a path connecting two TRIM in the same BZ, the 
Fermi energy inside the bulk gap will intersect these singly degenerate surface states either 
an even or odd number of times. When there are an even number of surface state crossings, 
the surface states are topologically trivial because weak disorder (as may arise through 
alloying) or correlations can remove pairs of such crossings by pushing the surface bands 
entirely above or below Ep. When there are an odd number of crossings, however, at least 
one surface state must remain gapless, which makes it non-trivial [3, S, Is]. The existence 
of such topologically non-trivial surface states can be theoretically predicted on the basis 
of the bulk band structure only, using the Z2 invariant that is related to the quantum 
Hall Chern number |24]. Materials with band structures with Z2 = +1 are ordinary Bloch 
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band insulators that are topologically equivalent to the filled shell atomic insulator, and are 
predicted to exhibit an even number (including zero) of surface state crossings. Materials 
with bulk band structures with Z2 = — 1 on the other hand, which are expected to exist in 
rare systems with strong spin-orbit coupling acting as an internal quantizing magnetic field 
on the electron system |[6j], and inverted bands at an odd number of high symmetry points 
in their bulk 3D BZs, are predicted to exhibit an odd number of surface state crossings, 
precluding their adiabatic continuation to the atomic insulator |10|, lul, I12I, ll3|. 

Such "topological quantum Hall metals" o], [lo|, 11 1 cannot be realized in a purely 2D electron 
gas system such as the one realized at the interface of GaAs/GaAlAs systems. 

In our experimental case, namely the (111) surface of Bio.gSbo.i, the four TRIM are 
located at T and three M points that are rotated by 60° relative to one another. Owing to 
the three-fold crystal symmetry (A7 bulk structure) and the observed mirror symmetry of 
the surface Fermi surface across = (Fig.2), these three M points are equivalent (and we 
henceforth refer to them as a single point, M). The mirror symmetry [E{ky) = E{—ky)] is 
also expected from time reversal invariance exhibited by the system. The complete details of 
the surface state dispersion observed in our experiments along a path connecting T and M 
are shown in Fig. 3a; finding this information is made possible by our experimental separation 
of surface states from bulk states. As for bismuth (Bi), two surface bands emerge from the 



bu 



k band continuum near F to form a central electron pocket and an adjacent hole lobe 



25I . l26l . l27l |. It has been established that these two bands result from the spin-splitting of a 



surface state and are thus singly degenerate [271, l28| . On the other hand, the surface band 
that crosses Ep at —k^ ~ 0.5 A~^, and forms the narrow electron pocket around M, is 
clearly doubly degenerate, as far as we can determine within our experimental resolution. 
This is indicated by its splitting below Ep between —k^ ~ 0.55 A~^ and M, as well as 
the fact that this splitting goes to zero at M in accordance with Kramers theorem. In 
semimetallic single crystal bismuth, only a single surface band is observed to form the 
electron pocket around M 0, [3^. Moreover, this surface state overlaps, hence becomes 
degenerate with, the bulk conduction band at L (L projects to the surface M point) owing 
to the semimetallic character of Bi (Fig.3b). In Bio.gSbo.i on the other hand, the states near 
M fall completely inside the bulk energy gap preserving their purely surface character at 



M (Fig.3a). The sur 



(unlike in Bi [25 



26 



'ace Kramers doublet point can thus be defined in the bulk insulator 



271, 



28j, 



29 



30|) and is experimentally located in Bio.gSbo.i samples to 
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lie approximately 15 ± 5 meV below Ep at k = M (Fig. 3a). For the precise location of this 
Kramers point, it is important to demonstrate that our alignment is strictly along the T — M 
line. To do so, we contrast high resolution ARPES measurements taken along the T — M line 
with those that are slightly offset from it (Fig.3e). Figs 3f-i show that with ky offset from the 
Kramers point at M by less than 0.02 A~^, the degeneracy is lifted and only one band crosses 
Ep to form part of the bow-shaped electron distribution (Fig. 3d). Our finding of five surface 
state crossings (an odd rather than an even number) between T and M (Fig.3a), confirmed 
by our observation of the Kramers degenerate point at the TRIM, indicates that these gapless 
surface states are topologically non-trivial. This corroborates our bulk electronic structure 
result that Bio.gSbo.i is in the insulating band-inverted {Z2 = —1) regime (Fig. 3c), which 
contains an odd number of bulk (gapped) Dirac points, and is topologically analogous to an 
integer quantum spin Hall insulator. 

Our experimental results taken collectively strongly suggest that Bio.gSbo.i is quite dis- 
tinct from graphene [isl, 3] and represents a novel state of quantum matter: a strongly 
spin-orbit coupled insulator with an odd number of Dirac points with a negative Z2 topo- 
logical Hall phase, which realizes the "parity anomaly without Fermion doubling". Our 
work further demonstrates a general methodology for possible future investigations of novel 
topological orders in exotic quantum matter. 

Note Added : In a very recent work we have successfully imaged the spin- 
polarization of the topological edge modes using high-resolution spin-resolved- 



ARPES 



3l|. 
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METHODS 



High resolution IPEM-ARPES data have been taken at Beamhnes 12.0.1 and 10.0.1 of the 
Advanced Light Source in Lawrence Berkeley National Laboratory, as well as at PGM Beamline 
of the Synchrotron Radiation Center in Wisconsin, with photon energies from 17 to 55 eV and 
energy resolution from 9 to 40 meV and momentum (k-)rcsolution better than 1.5% of the surface 
Brillouin zone. Data were taken on high quality bulk single crystal Bii_j;Sba; at a temperature of 
15 K and chamber pressures better than 8 x 10~^^ torr. Throughout this paper, the bulk bands 
presented are from those measured in the third bulk Brillouin zone to ensure a high degree of 
signal-to-noise contrast, and the kz values are estimated using the standard free-electron final state 
approximation. 
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Figure 1. Dirac-like dispersion near the L point in the bulk Brillouin zone. Selected 
ARPES intensity maps of Bio.gSbo.i are shown along three ^-space cuts through the L point of 
the bulk 3D Brillouin zone (BZ). The presented data are taken in the third BZ with = 2.9 
A~"^ with a photon energy of 29 eV. The cuts arc along a, the ky direction, b, a direction rotated 
by approximately 10° from the ky direction, and c, the kx direction. Each cut shows a A-shaped 
bulk band whose tip lies below the Fermi level signalling a bulk gap. The surface states are 
denoted SS and are all identified in Fig.2 (for further identification via theoretical calculations 
see Supplementary Information), d. Momentum distribution curves (MDCs) corresponding to the 
intensity map in a. f. Log scale plot of the MDCs corresponding to the intensity map in c. The 
red lines are guides to the eye for the bulk features in the MDCs. e. Schematic of the bulk 3D 
BZ of Bii_3;Sba; and the 2D BZ of the projected (111) surface. The high symmetry points F, M 
and K of the surface BZ are labeled. Schematic evolution of bulk band energies as a function of x 
is shown. The L band inversion transition occurs at x ~ 0.04, where a 3D gapless Dirac point is 
realized, and the composition we study here (for which x = 0.1) is indicated by the green arrow. 
A more detailed phase diagram based on our experiments is shown in Fig.3c. 
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Figure 2. Dispersion along the cut in the kz-direction. Surface states are experimentally 
identified by studying their out-of-plane momentum dispersion through the systematic variation 
of incident photon energy, a, Energy distribution curves (EDCs) of Bio.gSbo.i with electrons at 
the Fermi level (Ep) maintained at a fixed in-plane momentum of {kx=0.8 A^^, ky=0.0 A~^) 
are obtained as a function of incident photon energy to identify states that exhibit no dispersion 
perpendicular to the (lll)-plane along the direction shown by the double-headed arrow labeled "3" 
in the inset (see Methods for detailed procedure). Selected EDC data sets with photon energies of 
28 eV to 32 eV in steps of 0.5 eV are shown for clarity. The non-energy dispersive {kz independent) 
peaks near Ep are the surface states (SS). b, ARPES intensity maps along cuts parallel to ky taken 
with electrons at Ep fixed at kx = 0.8 with respective photon energies of hf = 22 eV, 29 eV 
and 35 eV (for a conversion map from photon energy to kz see Supplementary Information). The 
faint A-shaped band at hu = 22 eV and hu = 35 eV shows some overlap with the bulk valence 
band at L {hu = 29 eV), suggesting that it is a resonant surface state degenerate with the bulk 
state in some limited k-range near Ep. The flat band of intensity centered about —2 eV in the hv 
= 22 eV scan originates from Bi 5d core level emission from second order light, c, Projection of the 
bulk BZ (black lines) onto the (111) surface BZ (green lines). Overlay (enlarged in inset) shows 
the high resolution Fermi surface (FS) of the metallic SS mode, which was obtained by integrating 
the ARPES intensity (taken with hu = 20 eV) from —15 mcV to 10 mcV relative to Ep. The six 
tear-drop shaped lobes of the surface FS close to F (center of BZ) show some intensity variation 
between them that is due to the relative orientation between the axes of the lobes and the axis 
of the detector slit. The six-fold symmetry was however confirmed by rotating the sample in the 
kx — ky plane. EDCs corresponding to the cuts A, B and C are also shown; these confirm the 
gapless character of the surface states in bulk insulating Bio.gSbo.i- 
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Figure 3. The topological gapless surface states in bulk insulating Bio.gSbo.i. a, 

The surface band dispersion second derivative image (SDI) of Bio.gSbo.i along f — M. The shaded 
white area shows the projection of the bulk bands based on ARPES data, as well as a rigid shift 
of the tight binding bands to sketch the unoccupied bands above the Fermi level. To maintain 
high momentum resolution, data were collected in two segments of momentum space, then the 
intensities were normalized using background level above the Fermi level. A non-intrinsic flat band 
of intensity near Ep generated by the SDI analysis was rejected to isolate the intrinsic dispersion. 
The Fermi crossings of the surface state are denoted by yellow circles, with the band near —kx ~ 0.5 
A^^ counted twice owing to double degeneracy. The red lines are guides to the eye. An in-plane 
rotation of the sample by 60° produced the same surface state dispersion. The EDCs along T — M 
are shown to the right. There are a total of five crossings from T — M which indicates that these 
surface states are topologically non-trivial. The number of surface state crossings in a material 
(with an odd number of Dirac points) is related to the topological Z2 invariant (see text), b. The 
resistivity curves of Bi and Bio.gSbo.i reflect the contrasting transport behaviours. The presented 
resistivity curve for pure bismuth has been multiplied by a factor of 80 for clarity, c, Schematic 
variation of bulk band energies of Bii-^-Sb^^ as a function of x (based on band calculations and on 



[7j, [iTI]). Bio.gSbo.i is a direct gap bulk Dirac point insulator well inside the inverted-band regime, 
and its surface forms a "topological metal" - the 2D analogue of the ID edge states in quantum spin 
Hall systems, d, ARPES intensity integrated within ±10 meV of Ep originating solely from the 
surface state crossings. The image was plotted by stacking along the negative direction a series 
of scans taken parallel to the ky direction, e. Outline of Bio.gSbo.i surface state ARPES intensity 
near Ep measured in d. White lines show scan directions "1" and "2". f. Surface band dispersion 
along direction "1" taken with hv = 28 eV and the corresponding EDCs (g). The surface Kramers 
degenerate point, critical in determining the topological Z2 class of a band insulator, is clearly 
seen at M, approximately 15 it 5 meV below Ep. (We note that the scans are taken along the 
negative kx direction, away from the bulk L point.) h. Surface band dispersion along direction 
"2" taken with hv = 28 eV and the corresponding EDCs (i). This scan no longer passes through 
the M-point, and the observation of two well separated bands indicates the absence of Kramers 
degeneracy as expected, which cross-checks the result in (a). 
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FIG. 2: Hsieh et al., NATURE 452, 970 (2008). 
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FIG. 3: Hsieh et al., NATURE 452, 970 (2008). 
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